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PREFACE 


The  David  W.  Taylor  Lectures  were  conceived  to  honor  our  founder  in 
recognition  of  his  many  contributions  to  naval  architecture  and  naval 
hydrodynamics.  Admiral  Taylor  was  a  pioneer  in  the  use  of  hydrodynamic 
theory  and  mathematics  for  the  solution  of  naval  problems.  He  established 
a  tradition  of  applied  scientific  research  at  the  "Model  Basin"  which 
has  been  carefully  nurtured  through  the  decades  and  which  we  treasure 
and  maintain  today.  It  is  in  this  spirit  that  we  have  invited  Prof. 

Karl  Wieghardt  to  be  a  David  W.  Taylor  Lecturer. 

Prof.  Wieghardt  was  born  in  Vienna,  Austria,  in  1913.  From  1932 
to  1935  he  attended  the  Technical  High  School  in  Dresden,  Germany.  He 
began  studies  at  the  University  of  Gottingen  in  1935  and  received  the 
degree  of  Dr.  rer.  nat.  in  1938  for  studies  of  the  lift  distribution  on 
a  rectangular  airfoil.  In  1945  he  received  a  Dr.  habil  degree  for 
research  on  the  energy  equation  of  boundary  layers  and  earned  the  status 
of  university  lecturer.  From  1938  to  1949  he  was  employed  by  the 
Kaiser-Wilhelm  Institute  for  Flow  Research  (now  the  Max-Planck  Institute) , 
where  he  did  aerodynamics  and  general  flow  research.  In  1949  he  went 
to  Teddington,  England,  as  a  consultant  on  hydrodynamics  for  the  Admiralty 
Research  Laboratory.  Returning  to  Germany  three  years  later,  he  joined  the 
faculty  of  the  Institute  for  Shipbuilding  at  the  University  of  Hamburg, 
where  he  has  been  a  full  professor  since  1960.  He  is  a  member  and  past 
president  of  the  German  Society  for  Applied  Mathematics  and  Mechanics 
(GAMM)  and  a  past  member  of  the  Resistance  Committee  of  the  International 
Towing  Tank  Conference  (1957-1969). 


KINEMATICS  OF  SHIP  WAKE  FLOW 


INTRODUCTION 

.  For  the  usual  merchant  ship  at  constant  speed  in  calm  water  one  can  calculate 
the  wave  drag  approximately  and  even  the  three-dimensional  boundary  layer  up  to  say 
90%  of  the  ship  length.  But  we  do  not  know  where  and  how  the  boundary  layer  sepa¬ 
rates  and  mixes  with  the  general  secondary  flow  behind  the  ship  to  form  the  near¬ 
wake  flow  which  the  propeller  has  to  face  and  which  the  propeller  itself  influences. 
Isotachs  of  the  nominal  or  even  of  the  effective  wake  can  tell  only  part  of  the 
truth. 

Along  the  ship,  vorticity  is  produced  in  the  boundary  layer,  and  the  vector 
oo  =  rot  v  'lies  in  the  vertical  plane,  parallel  to  the  frame  and  perpendicular  to  the 
main  flow.  So,  a  vorticity  line,  i.e.,  a  line  parallel  to  the  vorticity  vector  at 
any  point,  would  just  go  round  the  hull  girthwise  at  the  parallel  midship.  Yet, 
further  on  in  the  wake,  one  usually  finds  what  looks  like  a  strong  longitudinal 
vortex  pair,  even  when  there  are  no  bilge  vortices  as  with  our  model.  Hence,  our 
starting  question  was:  how  are  these  vorticity  lines  bent  and  bundled?  Usual 
boundary-layer  tests  normal  to  the  hull  cannot  answer  this  question.  What  we  wanted 

was  a  general  survey  of  the  details  of  the  mean  flow  near  the  stern^ 

\ 

For  this  we  used  a  double  model  of  a  ship  with  block  coefficient  85  in  a  wind 

tunnel  (scale  1:95):  length  2.74  m,  breadth  0.404  m,  depth  0.148  m,  and  air  speed 

27  m/s,  giving  a  Reynolds  number  of  5  million.  The  model  was  hung  up  on  wires  in  a 

slotted  wall  test  section  (diameter  1.2  m,  length  6  m) .  The  pressure  distribution 

on  the  model  corresponded  very  well  to  that  of  potential  flow  as  calculated  by  the 

1* 

Hess-Smith  method,  except  for  the  stern  itself,  of  course. 

Measurements  were  made  in  12  vertical  planes,  x  =  const,  at  a  narrow  3-  by  3-mm 
grid  in  the  y-  and  z-directions  (sideward  and  upward).  By  linear  numerical  differen¬ 
tiation  in  the  y-  and  z-directions  and  by  graphical  differentiation  in  the  longi¬ 
tudinal  x-direction  all  nine  velocity  gradients  could  be  determined.  Figure  1  shows 
the  test  planes  and  their  sections.  In  some  planes  there  were  more  than  1,000  test 
points.  Basic  data  at  each  point  are  three  velocities  and  the  static  pressure, 
giving  nine  velocity  gradients  and  three  vorticity  components,  i.e.,  16,000  data 
in  one  single  plane. 


*A  complete  listing  of  references  is  given  on  page  63. 
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Obviously,  such  tests  can  be  made  only  by  on-line  measuring  with  a  computer. 

All  these  tests  were  made  by  Dr.  J.  Kux  who  also  had  introduced  Laser-Doppler- 
Velocimetry  (LDV)  to  the  Hamburg  Institute.  However,  the  test  probe  here  had  to  be 
as  simple  as  possible.  Hence,  a  calibrated  five-hole  tube  (diameter  3  mm)  was 
always  held  in  the  x-direction.  Only  tests  with  a  positive  pressure  in  the  front 
hole  were  evaluated.  So  we  have  no  data  not  only  next  to  the  body  but  also  in  the 
most  interesting  regions  with  the  strongest  cross  flow,  turbulence,  and  shear.  Yet, 
at  least,  now  we  know  where  hot  wire  or  LDV  tests  are  mostly  needed  as  a  supplement. 

It  is,  of  course,  impossible  to  state  precise  error  limits  for  all  test  data. 
Yet,  the  continuity  equation  renders  a  certain  control,  although  one  must  not  expect 
too  much  of  the  sum  of  three  differential  quotients  derived  from  measured  curves. 
Since  div  v  has  the  dimension  1/time,  we  divided  it  by  the  deformation  rate  d  as 
shown  in  Figure  2  vs  y  for  z  =  -1.5  mm  (z=0  is  the  horizontal  plane  through  the  axis 
of  the  nonexisting  propeller).  At  least  for  the  three  stern  planes  x  =  -38,  -19, 
and  0  mm  we  get  rather  similar  "curves."  That  is,  chance  errors  of  the  differenti¬ 
ation  must  be  small  compared  with  systematic  errors  of  the  tube  due  to  turbulence 
and  shear  at  comparable  side  distances  y.  For  example,  for  the  plane  x  =  0,  we 
know  from  LDV  tests  that  there  the  local  turbulence  degree  can  be  over  100%,  so  that 
the  rms  value  of  the  longitudinal  fluctuation  is  greater  than  the  mean  velocity. 
Besides,  any  tube  of  finite  size  straightens  the  flow  in  its  neighborhood.  Hence, 
in  particular  in  the  most  interesting  regions,  we  will  have  systematic  errors. 
However,  results  found  with  little  scatter  everywhere  in  the  flow  field  may  be 
trusted,  at  least  qualitatively. 

EXPERIMENTAL  FINDINGS 

Figures  2  through  11  show  the  secondary  flow  in  the  test  planes.  It  is  directed 
upward  and  inward  to  the  midship  plane  to  fill  the  hole  in  the  fluid  which  the  mid¬ 
ship  has  produced.  The  horizontal  plane  through  the  propeller  axis  is  at  z  =  0,  the 
keel  at  z  =  -44.5  mm.  All  velocities  are  normalized  by  U^,  the  velocity  of  the 
oncoming  flow.  At  the  end  of  the  parallel  midship,  at  x  =  -858  mm  in  Figure  2,  no 
bilge  vortices  have  formed.  The  blank  region  at  z  £  0  in  Figures  5  through  7 
indicates  a  separated  flow  region.  In  Figure  11  a  holding  strut  (diam  30  ram)  is  to 
be  seen. 
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It  is  amazing  how  smooth  and  repeatable  these  data  are.  Measurements  were  made 
line  after  line,  as  in  TV,  and  the  test  time  for  each  plane  was  more  than  a  week. 
Yet,  only  in  a  few  cases  did  differentiated  quantities  of  one  line  not  fit  those  of 
the  neighboring  lines.  Usually,  even  in  the  wake,  the  mean  flow  turned  out  to  be 
more  stable  and  stationary  than  could  be  expected. 

Isotachs  of  the  longitudinal  velocity  component  u  are  given  for  some  planes 
in  Figures  12  through  16.  The  distortions  of  a  few  lines  near  the  keel  are  caused 
by  the  wake  of  a  span  wire  under  the  keel  of  the  model  with  a  diameter  of  1.2  mm  at 
x  =  -380  mm.  Separation  seems  again  to  begin  near  z  =  0  (Figures  13  through  16). 

At  the  edge  of  this  zone  strong  velocity  gradients  are  to  be  expected. 

Figures  3  through  11  of  the  secondary  flow  suggest  the  formation  of  a  longi¬ 
tudinal  vorticity  component  =  w^  -  v Numerically  it  could  easily  be  found  by 
integrating  the  circulation  around  each  little  square  3  by  3  mm.  At  the  end  of  the 
parallel  midship  (Figure  2)  it  is  less  than  5/m  in  our  test  region,  i.e.,  at  least 
a  few  millimeters  off  the  model  wall.  (Since  we  have  divided  all  velocities  by  U 

J  CO 

the  dimension  of  our  oj  is  1/m  instead  of  1/s.)  Yet,  already  at  x  =  -157  mm  (Figure 
21)  the  maximum  of  is  about  20/m,  and  this  remains  almost  the  same  down  into  the 
wake  at  x  =  100  mm  (Figures  22,  23,  and  24). 

In  Figures  22  and  23  there  is  also  a  longitudinal  vorticity  In  the  opposite, 
negative  sense  up  to  -5/m.  To  explain  this,  one  may  imagine  the  double  model 
turned  by  90°  about  the  longitudinal  axis;  see  the  double  model  of  a  ship  with  a 

small  breadth  and  a  rising  keel  line  toward  the  stern.  Near  the  new  keel  (the 

former  waterline)  the  cross  flow  would  be  smaller  and  it  would  produce  negative 
vorticity  since  the  tests  were  now  made  on  the  port  instead  of  the  starboard  side. 

Since  the  longitudinal  vorticity  changes  so  little  downstream,  one  might, 
perhaps,  describe  it  by  tubes  of  potential  vortices.  In  the  wake  these  starboard 
and  port  vortices  attract  each  other  and  the  originally  V-shaped  vorticity  region  is 
folding  up. 

But  this  would  not  be  the  main  point.  For  it  turns  out  that  the  cross  compo¬ 
nents  of  vorticity  and  U)^  are  mostly  much  larger  than  the  longitudinal  component 
oj^.  In  Figure  25  at  x  =  -73  mm  the  magnitude  of  the  complete  vorticity  vector  /w/ 
has  a  maximum  value  of  about  80/m  in  the  test  region.  At  x  =  -19  mm  (Figure  26)  it 

is  still  80/m  but  now  at.  the  edge  of  the  separated  region  off  the  wall.  Further 

downstream,  in  the  wake,  the  maximum  has  decreased  to  less  than  50/m  (Figure  27). 
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Incidentally,  such  a  region  with  high  vorticity  (and  Reynolds  stresses)  must 
not  be  confused  with  a  strong  potential  vortex  with  high  underpressure  in  its  core. 
Here,  the  static  pressure  field  is  rather  even  as  Figure  28  shows  for  x  =  -19  mm. 
Also  in  the  wake,  at  x  =  50  mm  (Figure  29) ,  the  pressure  changes  only  little  in 
spite  of  the  high  cross  vorticity  of  the  mean  velocity  field. 

Another  way  to  become  acquainted  with  this  intricate  three-dimensional  mean 
flow  is  to  ask  for  the  angle  between  the  velocity  and  vorticity  vectors.  Since 
usually  the  velocity  is  directed  mainly  in  the  longitudinal  x-direction,  this  angle 
is  also  approximately  the  one  between  x  and  to.  For  x  =  -19  mm,  Figure  30  shows  this 

angle  ^  v,  to.  An  acute  angle  is  found  only  in  the  hollow  groove  of  the  section 
with  the  maximum  of  longitudinal  vorticity  to^.  Only  there  is  a  kind  of  longitudinal 
vortex  roll.  Everywhere  else,  especially  in  the  region  with  /to/^  ,  velocity  and 
vorticity  are  almost  perpendicular  to  each  other.  This  is  also  true  for  the  wake 
at  x  =  50  mm  (Figure  31)  where  the  angle  between  v  and  to  becomes  small  only  near  the 
vertical  midship  plane  due  to  symmetry  reasons. 

That  velocity  and  vorticity  are  mostly  perpendicular  to  each  other,  not  only  in 
any  two-dimensional  flow  or  in  a  boundary  but  also  after  separation  in  the  wake  and 
at  greater  distances  from  the  wall,  becomes  quite  trivial  when  looking  at  the  time 
average  of  the  Navier-Stokes  equation 


a  =  grad  v  / 2 


v  x  to  *  -  —  grad  p  -  V  rot  to  +  K(Re.str.) 


where  K(Re.str.)  represents  the  resultant  force  due  to  Reynolds  stresses.  Contours 

of  /a/  at  x  =  -19  are  shown  in  Figure  32.  In  wake  flow,  because  there  are  still_ 

->2 

velocity  differences  in  particular  in  the  longitudinal  component,  the  term  grad  v  /2 
is  still  existent.  On  the  other  hand,  the  static  pressure  field  is  almost  level  and 
the  viscous  term  for  the  mean  flow  is  negligible.  Farther  on,  far  away  from  the 
hull,  Reynolds  stresses  and  their  resultant  forces  are  also  small.  Hence,  the 
second  term  v  x  u  must  almost  cancel  the  first  one,  and  this  implies  that  the  angle 
between  v  and  to  is  not  too  small. 

When  one  does  trust  in  all  these  data,  one  may  add  the  pressure  gradients  to 
the  acceleration  components  to  get  the  resultant  forces  of  the  Reynolds  stresses. 

An  example  is  given  in  Figure  33:  the  longitudinal  force  component  in  the 
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at  various 


horizontal  plane  z  =  -1.5  mm  vs  the  distance  from  the  hull  y  -  y^  ^ 
planes  of  constant  x.  At  x  =  -63  mm  the  boundary  layer  is  still  adjacent  as  the 
turbulent  friction  force  is  adverse  to  motion,  i.e.,  in  the  negative  x-direction. 

But  at  x  =  -53  mm  this  force  has  disappeared  or  is  even  in  the  positive  x-direction 
near  the  wall.  Hence,  in  this  z-plane,  separation  seems  to  occur  between  x  =  -63 
and  -53  mm.  This  would  be  difficult  to  find  out  from  velocity  profiles  alone  as 
given  in  Figure  34:  u  vs  wall  distance  Ay  =  y  -  y^^ >  which  is,  of  course,  not  the 
correct  normal  to  the  hull  (the  static  pressure  is  fairly  constant  in  each  x-plane) . 
Flow  directions  projected  to  the  horizontal  plane  z  =  -1.5  mm  are  shown  in  Figure 
35. 

Another  example  of  Reynolds  forces  is  given  in  Figure  36  where  lines  of  =  -2, 
-4,  and  -8/m  are  shown  in  the  plane  x  =  -19  mm.  In  some  points  near  the  keel  all 
three  force  components  are  indicated.  Near  the  keel  the  Reynolds  forces  are  directed 
toward  the  keel  to  compensate  for  the  opposite  gradients  there  (cf  Figure  28)  and 
the  inertial  forces.  In  principle,  one  could  even  try  to  check  the  hypothesis  of 
eddy  viscosity  locally;  but  then,  the  velocity  components  would  have  to  be  dif¬ 
ferentiated  twice. 

At  last,  a  remark  on  a  similarity  rule.  Professor  Mori,  Hiroshima,  has  told  us 
that  he  had  also  found  high  cross  vorticity  near  the  stern  of  a  towed  3-m  model  as 
published  in  1975.  Some  of  his  results  are  plotted  in  Figure  37  where  /w*/  =  L/2 
j oo | / with  L  =  model  length.  To  compare  these  tests  in  water  with  ours  in  air,  we 
use  U  and  the  boundary-layer  thickness  6  at  the  end  of  the  model  corresponding  to 
the  end  of  a  flat  plate  with  the  same  length 

6  =  0.085  L  Rn-0,1  +  3%  for  106 * *  <  Rn  £  109 

3  |  -*  | 

as  derived  from  measurements  by  Winter  and  Gaudet.  Then  the  maximum  for  oj  6/U  is 
4  for  our  double  model  in  air  and  2.6  for  the  Japanese  model  in  water  with  a  coarser 
test  grid.  Hence,  at  least  for  the  order  of  magnitude  one  may  expect  for  a  ship 

9 

with  L  =  200  m,  U  =  7.1  m/s,  Fn  =  0.16,  Rn  =  1.2  •  10  ,  and  6  =  2.1  m  a  cross 

4 

vorticity  up  to  10/s. 
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* 
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TWO  VORTICITY  NUMBERS 

The  most  spectacular  test  result  is  the  region  with  high  cross  vorticity.  Yet, 
when  can  we  speak  of  high  or  low  vorticity?  Strictly  speaking  this  question  does  not 
make  sense  because  vorticity  has  a  dimension,  namely  1/s  or  here,  after  division  with 
U  1/m.  Only  with  solid  bodies,  quantities  such  as  strain,  torsion,  or  shear  are 
dimensionless;  with  fluids  they  are  rates.  Hence,  C.  Truesdell^  introduced  a  local, 
kinematic  vorticity  number  where 

W  =  |  to  |  /  d  =  magnitude  of  vorticity  over  deformation  rate,  (d^=co^+2  div  a), 

K 

— ► 

a  =  acceleration  in  stationary  flow,  or 

-+2  .  “V 

a  =  grad  v  /2  -  v  x  or  written  in  full 


<  ■ 


d  =  2(u  +v  +w  )  +  (w  +v  )  +  (u  +w  )  +  (v  +u  ) 

xyz  y  z  zx  xy 

1  V  J 

\  — -  M  1  ■  — - — f 


strains 


"shears" 


For  example,  for  a  straight  vortex  line  with  peripheral  velocity  proportional  to 
(radius)0  this  number  becomes  =  |n+l/n-lj .  Hence,  for  the  potential  vortex  with 
n  =  -1,  =  0,  of  course,  but  for  a  rotating  rigid  body  without  any  deformation 

n  =  1  and  W  -»■  °°.  In  pure  shear  flow  and  hence  in  any  tube  flow  (laminar  or  turbu- 

-*• 

lent)  W^  =  1  exactly  because  there  is  no  acceleration  and  no  div  a.  In  a  boundary 
layer  is  practically  also  equal  to  1,  but  only  approximately. 

Applying  this  measure  to  the  mean  velocity  field  of  our  turbulent  stern  flow 
yields  the  lines  W  in  Figures  38,  26,  and  39  for  the  planes  x  =  -73,  -19,  and 
50  mm.  The  region  with  the  highest  |w] -values  has  become  "invisible,”  since  here 

W  is  everywhere  between  0.75  and  1;  i.e.,  almost  all  deformation  is  due  to  vorticity 

K  ^  -+2 
and  div  a  is  positive  but  small  compared  with  w  . 

On  the  other  hand,  in  the  hollow  groove  of  the  sections  a  longitudinal  vortex 

roll  is  formed  and  marked  by  a  maximal  value  of  Wv  -  2  at  x  *  -19  (Figure  26). 

^  -*• 

Here,  there  is  more  rotation  than  deformation  of  the  fluid  particles  and  div  a  is 
negative. 

At  greater  wall  distances  and  in  particular  later  in  the  wake,  W^,  decreases  to 
and  below  0.5  at  the  edge  of  the  test  range.  Obviously,  the  fluid  particles 


l 


originally  contaminated  with  vorticity  produced  in  the  boundary  layer  have  been 
diluted  by  the  nonrotating  particles  swept  into  the  wake  from  outside  by  the 
secondary  transverse  flow. 

C.  Truesdell^  was  also  looking  for  a  simple  local  measure  indicating  the  magni¬ 
tude  of  error  one  would  make  in  calculating  a  laminar  flow  by  neglecting  vorticity 
or  assuming  potential  flow.  As  a  dynamic  vorticity  number  he  defined 

I  VX(jj| 

w° "  -2/2i 

I  +  grad  v  /  2  j 

or  for  stationary  flow  simply 

WQ  =  jvxwj/|grad  v2/2| 

For  a  straight  vortex  line  with  peripheral  velocity  proportional  to  r°,  W^  = 
|(n+l)/n|,  i.e.,  for  the  potential  vortex  with  n  =  -1,  WQ  =  0  for  r  >  0,  yet  for  a 
fluid  rotating  like  a  rigid  body,  n  =  1  and  W  =  2.  (In  incompressible  airfoil  or 
propeller  theory  the  free  vortices  leave  along  streamlines  so  that  v||w  and 
Wp=0  everywhere.)  Where  there  is  no  acceleration,  e.g.,  in  stationary  tube  flow, 

WD  *  1.  Again,  in  a  boundary  layer  W^  is  nearly  1. 

The  distribution  of  W^  in  our  stern  flow  is  shown  in  Figures  40,  41,  and  42 
for  x  =  -73,  -19  and  50  mm.  In  the  region  with  the  highest  |co| -values,  WQ  is  also 
about  1,  similar  to  W^,.  Again,  only  in  the  longitudinal  vortex  roll  does  reach  a 
maximum  of  up  to  2  (Figure  42).  Yet,  it  should  be  mentioned  that  the  calculation  of 
Wp  showed  more  scatter  than  that  of  Wj,.  Hence,  quantitatively  these  maximal  values 
of  W^  are  probably  not  as  reliable  as  those  of  Wv. 

D  ix 

However,  it  may  be  concluded  that  in  a  wide  region  W^  is  near  to  1  and  at  the 
edge  of  stern  flow  and  wake,  W^  also  decreases  at  least  to  about  0.5  in  the  test 
region,  as  does  W^. 

TRANSITION  FROM  ROTATIONAL  TO  IRROTATIONAL  FLOW 

The  decay  of  W^  and  W^  at  greater  distances  from  the  body  seems  quite  trivial. 
Inside  a  boundary  layer  or  wake  the  deformation  consists  mainly  of  rotation,  but  in 


7 


r 

^  - 


(the  potential  flow  outside  there  is  only  deformation  (and  change  of  kinematic  energy 
-*2 

"  pv  /2)  but  no  longer  rotation.  Yet,  surprisingly,  in  this  respect  the  test  results 

are  more  precise  than  all  boundary-layer  calculations,  simply  because  the  usual 
1  assumptions  of  boundary-layer  theory  (e.g.,  9v/9x<<3u/9y)  do  not  hold  any  longer  at 

the  edge  of  the  boundary  layer.  There,  the  asymptotic  transition  into  the  outer 
potential  flow  is  described  correctly  only  to  the  first  order  for  the  velocities  but 
not  for  the  velocity  gradients.  Already  Blasius'  solution  for  the  flat  plate  gives 
div  a  £  0  everywhere  and  hence  £  1  for  any  wall  distance. 

The  same  deficiency  is  found  with  solutions  for  free  plane  flows  of  a  boundary 
layer  nature,  meaning  that  the  region  of  space  in  which  a  solution  is  being  sought 
does  not  extend  far  in  a  traverse  direction  as  compared  with  the  main  direction  of 
i  flow,  whether  these  flows  are  laminar  or  turbulent. 

such  stationary  plane  flows  with  similarity  solutions  the  stream  function 

ip  =  C  xnF(n),  resp  Uy  -  C  xRF(n)  for  wake  flow,  n  =  ayx^ 
divergence  of  acceleration  is 

2  div  a  =  4  C2o2x2 (n+k_1)  .  { (n+k)2F,2+k(k+l)nF'F"-n(n-l)F*F"} 


For 


is 


6,7 


Then  the 


and  W  = 


|/d  =  l//l+2(div  a)/w2 


\ 


)  ■ 

.. 


Examples  are  as  follows: 

1.  Boundary  layer  along 
2  div  a  =  U2x_2(F-nF')  F"  £  0 

2.  The  smoothing  out  of 
turbulent  jet:  n  *  1,  k  =  -1 
function  F(h)  or  any  velocity 

3.  Laminar  free  jet:  n 
2  F  F"  <  0  and  WR  >  1  for  n  > 


a  flat  plate:  n  =  1/2,  k  =  -1/2, 
and  £  1  everywhere. 

a  velocity  discontinuity  or  the  boundary  of  a  plane 
and  div  a  =  0  and  W  =  1  everywhere  and  even  for  any 
profile  F' . 

=  1/3,  k  =  -2/3,  F  =  tanh  n,  2  div  a  ~  F'2  -  2r)F'F"  + 
0.794. 
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4.  Turbulent  free  jet:  n  =  1/2,  k  =  -1;  approximation  functions 
tanh  n  gives  W  >  1  for  n  >  0.658 

K  — 


_i 

/e 


dt  W  >  1  for  n  >  0.571, 

K  —  — 


as  shown  in  Figure  43. 


for  ip: 


5.  Plane  wake,  laminar  or  turbulent: 


n  =  0, 


k  =  -1/2  and 


dn. 


Similarity  is  found  experimentally  for  Ux/4e  >  600  with  £  =  kinematic  or  eddy 

viscosity;  then  2(div  a)/ui^  <  0.001  or  >  0.999,. 

K  j 

Unfortunately,  for  the  transition  from  rotational  to  irrotational  flow  there  is 

only  one  exact  example  known,  namely  the  stagnation  point  flow.  It  is  again  a 

similarity  solution  and  for  plane  laminar  flow  n=l,  k=0,  ip  =  V  £F(n)  with 

5  =  x/a/v,  n  =  y/a/v  and  F'"+F4F"-F' ^+1  =  0.  Here,  2  div  a  =  4  a^F'^  >  0  or  W„  £  1. 

As  shown  in  Figures  44  and  45,  W  and  W_  are  both  1  at  the  wall,  and  they  decrease 

K.  D 

monotonously  with  wall  distance  to  zero.  The  axisymmetr ic  stagnation  point  flow 
also  shows  this  expected  behavior. 

Another  exact  solution  of  which  one  might  think  in  this  connection  is  the  one 

for  the  flow  in  a  convergent  or  divergent  channel.  There  is  only  a  radial  velocity 

2  ->-2  2  4 

v  f(a)/r  (in  polar  coordinates)  with  oj  =  -v  f'(a)/r  and  2  div  a  =  4  v  f  ( a)/r  >.  0. 

At  the  walls  f  =  0,  f '  ?  0,  WR  =  1;  at  the  center  plane  f  ^  0,  f '  *  0,  WR  =  0. 

(Only  in  the  limiting  case  of  flow  between  two  parallel  walls  is  there  pure  shear 

flow  without  acceleration  and  W  =  1  everywhere.)  Yet,  there  is  no  transition  to 

irrotational  flow.  At  the  center  plane  Wv  must  vanish  for  symmetry  reasons,  since  w 

is  changing  sign  there,  just  as  it  does  at  the  center  plane  of  a  jet  or  wake. 

A  more  pertinent  example  is  the  round  laminar  jet:  in  spherical  polar  coordi- 

nates  ip  =  vR*f(n=  cos  0)  with  f  =  2(l-n  )/(e+l-n)  and  e  «  1.  Figure  46  gives  the 

velocity  distribution  in  usual  cylindrical  coordinates,  i.e.,  the  axial  component 

u  =  vD  cos  9  -  v0  sin  0  vs  radial  distance  from  the  axis  r  =  R  tan  0  over  b  =  R 

tan  0.  with  u/u  =  0.5  at  r  =  b.  The  arbitrary  constant  £  is  assumed  to  be  £  =  0.01 
0  2 

as  in  Reference  7,  corresponding  to  M/pV  =  3282.5  with  M  =  momentum  of  the  jet; 
then  9*  =  5.184°.  is  about  1  in  a  wide  range,  viz.,  WR  =  1+5%  for  0.3  <  r/b  < 
4.6  where  u/u  falls  from  0.93  to  0.01.  Only  at  greater  axial  distances  W  decays; 

O 

e.g.,  W  -  0.5  for  r/b  =  6.67  where  u/u  =  0.0037.  In  this  special  case  the  approxi- 

K  O 

mate  "boundary  layer  solution"  corresponds  precisely  to  the  limit  of  the  exact  so- 

2  2 

lution  for  £  -*•  0.  For  small  angles  0  near  the  axis  f  becomes  £  / (1+5  / 4)  with 
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£  =  0/2/e  =  0/4v  (3/tt  M/p)1/2,  since  for  e  -*■  0  M/pv2  -*•  32  ir/3e.  This  yields  u/u  = 

2  2  ° 
l/[l+(/2-l)(r/b)  ]  .  Within  the  accuracy  of  drawing,  both  curves  for  u/u  and  W 

o  K 

from  this  solution  are  indiscernible  from  those  of  the  exact  solution,  though,  for 

example,  at  r/b  =  6.67,  the  approximate  solution  gives  again  W„  =  0.3  but  u/u  = 

K.  O 

0.0026  instead  of  0.0037. 

For  turbulent  flow  there  are  no  exact  solutions  and  here  the  transition  from 
rotational  to  irrotational  outer  flow  is  obscured  by  intermittency  anyway.  Farther 
on,  the  vorticity  numbers  W^,  and  derived  from  the  acceleration  a  and  rotation 
of  the  mean  velocity  field  are,  perhaps,  somewhat  artificial  physical  quantities. 
Possibly,  it  would  be  more  appropriate  to  use  mean  values  of  the  real  acceleration 


a*  and  rotation 


/+2  +.2 
y  to  +u)  ; 


with  v  and  to  velocity  and  vorticity  fluctuation 


grad  (v2+v'2  )  /2  -  (v+v')  x  (uH-to") 


->.2 


a  +  grad  v  / 2  -  v'xu'  = 


grad  p  -  v  rot  to 


Then  the  divergence  of  the  mean  of  the  real  acceleration  a*  gives,  as  for  laminar 
flow,  simply 

div  a*  =  -  Ap/p 

At  higher  Reynolds  numbers  the  viscous  term  V  rot  to  is  negligible,  as  has  been 
checked  in  our  tests;  then,  there  exists  an  acceleration  potential  for  a*  namely 
-  p/p  +  const.  However,  in  order  to  determine  "real"  vorticity  numbers,  e.g.,  WRA 


<  ■ 


W 


K* 


1/  V  l+2(dlv  a,)/  (w2-Ko 


-2) 


one  would  have  to  know  not  only  the  mean  velocity  field  in  all  details  but  the 
complete  tensor  of  Reynolds  stresses  as  well. 
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( 


1  +  3 


See  the  appendix  for  the  deformation  rate  in  potential  flow. 


SUMMARY 

The  stern  flow  and  the  near  wake  of  the  double  model  of  a  full  ship  in  a  wind 
tunnel  have  been  measured  at  so  many  points  that  not  only  the  mean  velocity  field 
but  also  all  nine  velocity  gradients  could  be  determined  with  sufficient  accuracy. 

Not  only  in  the  boundary  layer  still  adjacent  to  the  hull  but  also  after  separation 
and  in  the  near  wake  the  velocity  and  vorticity  vectors  are  mostly  nearly  perpendicu¬ 
lar  to  each  other.  There  is  a  region  with  especially  strong  cross  vorticity  of  the 
order  of  up  to  10/s  to  be  expected  on  a  large  ship. 

Two  vorticity  measures  of  C.  Truesdell,  W„  and  W  ,  are  used  to  analyze  this 

K  D 

turbulent  flow.  There  is  a  wide  region  with  W„  and  W  near  to  one,  and  outward  to 

the  edge  of  stern  flow  or  wake  both  vorticity  numbers  decay.  This  transition  from 

rotational  to  irrotational  flow  is  not  described  by  known  similarity  solutions  for 

laminar  or  turbulent  jets  and  wakes  because  the  usual  boundary  layer  assumptions  do 

not  hold  at  the  edge.  Only  the  exact  solution  for  laminar  stagnation  point  flow 

yields  correctly  WR  Q  -*•  0  for  increasing  wall  distance.  Since  the  region  with 

decreasing  W  is  very  large  for  near-wake  flow,  as  shown  by  these  tests,  an 
K ,  U 

analytical  study  of  such  transition  from  inner  rotational  to  outer  potential  flow 
seems  necessary  for  a  realistic  description  of  wake  flow. 


11 


■y 


20 


40 


60 


y  60 


X 

-38 

+ 

19 

• 

0 

O 

50 

& 

U*  *►  Vy  ♦  Wz 

2  (ux  t  vj -V  w|  )  +  ( Wy  4  vz)e 
+  {u,  ♦  W.)*  4  (v„  +  UyJ* 


Figure  lb  -  Deviations  of  Test  Data  from  Continuity  Equation 
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Figure  2  -  Cross  Flow  in  Plane  x  «  -858 
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Figure  6  -  Cross  Flow  in  Plane  x  =  -38 
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Figure  7  -  Cross  Flow  in  Plane  x  =  -19 
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Cross  Flow  In  Plane  x  =  0 
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Figure  9  -  Cross  Flow  in  Plane  x  =  50 
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Figure  11  -  Cross  Flow  in  Plane  x  =  200 
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Figure  14  -  Isotachs  of  Longitudinal  Velocity  u  in  Plane  x 
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x=-19.mm 


Figure  16  -  Isotachs  of  Longitudinal  Velocity  u  in  Plane  x  =  -19 
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tant  Angle  between  Velocity 


Figure  36  -  Force  Components  K  Resulting  from  Reynolds  Stresses:  Lines  for 
Kj,  -  const  and  Arrows  for  Cross  Component  +  in  the 
Vertical  Plane  z  *  -1.5  at  x  »  -19 
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Magnitude  of  Vorticity  |co*|  =  |w|^/U^  at  x/L  =  0.95  after  [2] 
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Stagnation  -  point  flow  ••  streamlines  ;  isobars,  vetocity  profile. 
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Laminar  Stagnation  Point  Flow:  Vorticity  Numbers  W„  and  W, 


Figure  46  -  Round  Laminar  Jet:  Velocity  Profile  and  vs  Distance  from  Axis 


APPENDIX 

REMARKS  ON  THE  DEFORMATION  RATE  IN  POTENTIAL  FLOW 

To  get  used  to  the  deformation  rate  in  incompressible  fluid,  it  might  be  help¬ 
ful  to  consider  simple  examples  from  plane  potential  flow,  although  here  WK  =  WD  = 
of  course.  In  general,  a  plane  flow  with  a  stream  function  ip  has  a  deformation 


rate  d: 


d^  =  +  2  div  a  =  +  4  •  (ip^  -ip  ip  ) 

T  xy  xx  yy 


In  a  potential  flow  with  w  =  -  Alp  =  0  and  F(z=x+iy)  =  +  iip 


d^  =  2  div  a  =  4(<J>^  -<p  <f>  ) 

xy  xx  yy 


=  4(t|T  ip  ) 
xy  xx  yy 


=  4  F"  F”  =  (2  |  F”  (  )  2  _>  0 

i.e.,  the  deformation  rate  is  simply  twice  the  magnitude  of  F" . 
This  gives  for  point  singularities  such  as 


source  F  =  -5—  In  z,  d  =  jql/irr^ 

2tt 


F  2 

vortex  F  =  In  z,  d  =  J  F )  / Tir 


dipole  F  =  m/(2Tiz),  d  =  |m  |  /  Trr  “ 


The  addition  of  a  parallel  flow  AF  =  U  z  does  not  alter  F" .  Hence,  the  lines  of 

constant  deformation  rate  are  circles  also  for  the  halfbody  or  the  circular  cyl- 

2 

inder.  For  inviscid  stagnation  point  flow  F  =  a/2  z  ,  the  deformation  rate  is  the 
same  everywhere,  d  *  2 | a [ . 
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After  conformal  mapping  £(z )  one  has 


with  R_  =  radius  of  curvature.  Again,  this  is  only  true  for  an  elliptical  strut. 

3 

Yet,  it  is  mentioned  here  because  K.  Oswatitsch  proposed  as  a  vorticity  number  in 
plane  flow  (including  compressible  flow)  the  product  of  local  vorticity  and  radius 
of  curvature  of  the  streamline  over  velocity:  |a)|  •  R^/U^.  For  an  elliptical  pro¬ 
file  this  is  exactly  proportional  to  W  =  M  /d  on  the  profile  itself. 

K. 
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Elliptical  Strut  in  Potential  Flow  at  Angle  of  Incii 
Streamlines  and  Lines  for  Constant  Deformation  Rate 
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